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Abstract: In this paper, we considered an inflationary model that effectively behaves 
as a modified Chaplygin gas in the context of quintessence cosmology. We reconstructed 
the inflaton potential bottom-up and using the recent observational data we hxed the free 
parameters of the model. We showed that the modified Chaplygin gas inspired model is 
suitable for both the early and the late time acceleration but has shortcomings between 
the two periods. 
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1 Introduction 

Chaplygin gas can be thought of as a perfect fluid with a bulk viscosity where the bulk 
viscosity is the sum of powers of density. Early solutions of this type of inflationary and 
deflationary models were studied in the context string theory [1, 2]. When the solution of 
bosonic d-brane in (d + 1,1) space-time has been studied classically [3] the name Chaplygin 
gas resurfaced in the context of cosmology. One can begin from Nambu-Goto action in the 
light-cone gauge and solve the momentum constraint to obtain the Chaplygin gas equation. 
Later, it was found that other forms of matter on the brane, such as Chaplygin gas, can 
be introduced [4], where a supersymmetric generalization is only possible for a fluid of 
the Chaplygin gas type [5]. Friedmann-Robertson-Walker (FRW) cosmological model with 
Chaplygin gas was proposed as an alternative to quintessence [6]. The Chaplygin gas 
equation of state, 

B , , 

P=--. (1.1) 

leads to a homogenous cosmology with, 

P{a) = + ( 1 - 2 ) 

where C is an integration constant. Therefore, it behaves like dark matter in the early uni¬ 
verse (a —>■ 0) and like dark energy in the late universe (a —>• oo). This realization led to a 
proposal [6] of inhomogeneous Chaplygin gas which can unify dark matter and dark energy 
[7]. The density perturbations and the mass power spectrum for a universe dominated by 
the Chaplygin gas have also been studied [8, 9]. The supernova type la observational data 
have been fitted using a model with the Chaplygin gas together CDM [10]. It has been 
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shown that the predicted age of the universe in the context of Chaplygin gas CDM models 
is in agreement with the most recent age estimates of globular clusters [11]. 


In order for the universe to shift from a dust dominated era to a de Sitter phase, the 
Chaplygin gas equation of state was generalized to [6] , 

P = -^, ( 1 - 3 ) 

with the following scale factor-dependent energy density [12], 

1 

l + o: 

( 1 - 4 ) 


p{a) = 


B + 


C 


a3(i+“) 


where C is an integration constant and a > 0 is a free parameter. This model is called the 
generalized Chaplygin gas model [13]. It is clear that a = 1 reproduces the pure Chaplygin 
gas model. Constraints on the generalized Chaplygin gas from supernovae observations 
obtained for the first time by the Ref. [14]. Also, statistical analysis of the 2dFGRS 
data yields to finding that very small (close to zero) and very large values (a S> 1) of 
the equation-of-state parameter a are preferred [15]. In the recent work [16] the non¬ 
linear regime of unified dark energy models, using the generalized Chaplygin gas cosmology 
has been studied. Recent observational constraint on generalized Chaplygin gas model 
indicated that B = 0.73^0 og, a = —O.OOIq}^ lu level and B = 0.73lQQg, a = —0.09^0 ig 
at 2(7 level [17] which agree with the claim of the Ref. [15] for the preferred small value 
of a while, interestingly, negative value of a was seen to be possible. Choosing negative a 
changes nature of the generalized Chaplygin gas to some barotropic matter with p = , 

where /3(= —a) > 0. In the other effort [18] using Markov Chain Monte Carlo method 
it is found that a = 0 . 001261 qqq)’ 26 %°ooi 26 * ^ ~ *^•’^'^^^0 0161^0 033 ^ While there are 

several attempt to fix B using observational data, some authors suggested B as variable 
parameter [19, 20] and the resulting model called varying generalized Chaplygin gas or new 
generalized Chaplygin gas model. There is also the possibility to consider viscosity in the 
generalized Chaplygin gas model [21]. 

The generalized Chaplygin gas equation of state can be extended to the following form 


[ 22 ], 


P = 


B 

l+CJ 


- 1 + +1 


B cj 

1+uj) 


(1.5) 


which is called the generalized cosmic Chaplygin gas with 0 > w > —I, and I being a 
positive definite constant which can take on values larger than unity [23]. There is also a 
class of equations of state interpolating between ordinary fluids with high energy densities 
and generalized Chaplygin gas fluids with low energy densities which is called the modified 
Chaplygin gas [24] . Soon after that viscosity was introduced to the modified Chaplygin gas 
[25]. Finally further extensions were introduced as modified cosmic Chaplygin gas [26-28], 
and extended Chaplygin gas [29-33]. The modified Chaplygin gas is subjects of this paper 
to study inflation. 

Inflationary scenarios including generalized Chaplygin gas have been considered already 
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[34-36]. Recently, the generalized Chaplygin gas inflation in the light of Planck [37] and 
BICEP2 [38] has been studied and has been shown that the generalized Chaplygin gas 
is not a suitable candidate for inflation in the context of Einstein theory of gravity [39]. 
Now, we would like to investigate the consequences of Inflation obeying the equation of 
state of the modified Chaplygin gas so that it can explain both the primordial and late 
time acceleration. We will use recent observational data to fix parameters of the modified 
Chaplygin gas model, in other words we would like to construct the inflaton potential in 
order to get a viable mechanism during both of the accelerated stages. 

The paper is organized as follows. In next section we introduce modified Chaplygin gas 
inflation and obtain the scale factor dependence of the scalar field and construct the inflaton 
potential by solving Eriedmann equations. Then, in section 3, we calculate the tensor to 
scalar ratio, slow-roll parameters, scalar spectral index by using the inflaton potential and 
use Planck data to hx the model parameters. Using hxed parameters we investigate some 
cosmological parameters in section 4. In section 5 we briefly examine the effects of varying 
bulk viscosity in modified Chaplygin gas. In section 6 we give our conclusions and suggest 
some future work. 


2 Modified Chaplygin gas infiation 

Although the Chaplygin gas models were introduced to explain late time acceleration with¬ 
out dark energy, the primordial acceleration which is believed to be driven by a scalar field, 
Inflaton, can also be described by the same modified Chaplygin gas model. Therefore we 
would like to describe both the primordial and recent inflationary phase with a scalar field, 
that we can still call inflaton, for which the equation of state is the one that the modified 
Chaplygin gas obeys. 


The MCG is dehned with the following equation of state [40], 


p = Ap 



( 2 . 1 ) 


where A, a, and B are free parameters of the model. The case of A = 0 recovers generalized 
Chaplygin gas EoS, and A = 0 together a = 1 recovers the original Chaplygin gas EoS. The 
special case of a = 0.5 has been studied in the Ref. [41]. Ref. [42] where it was concluded 
that the best fitted parameters are A = —0.085 and a = 1.724. Other observational con¬ 
straints on the modified Chaplygin gas model using Markov Chain Monte Carlo approach 
resulted into A = 0.00189 ^'^qQQ^ gg, a = 0.1079lQ2f39 at Icr level and A = 0.00189lQQQg®5, 
a = 0.1079^02911 at 2(T level [43]. Other interesting constraints may be also be found based 
on cosmic growth [44]. 


All of the observational constraints given above can be summarized with a single 
statement; parameter A may be a small positive or small negative number. Many of these 
analysis were done with a very strong simplification; Chaplygin gas being the only matter 
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in the universe. It is relatively safe to make that simplihcation when considering a non¬ 
zero parameter B] since for the final stages of the universe, with decreasing energy density, 
Chaplygin gas is going to dominate the universe and will determine the Hubble parameter. 

We would like to constrain the parameter A by assuming that the same scalar is caus¬ 
ing the primordial inflation unlike the studies above. Primordial inflation is described by 
quasi de Sitter geometry with the equation of state being w = —1 + e where e is a positive 
infinitesimally small number. Therefore, one can trivially infer from equation (2.1) that 
parameter A should be very close to -1 since that part of the pressure will dominate during 
the early universe. And the peculiar form of equation (2.1) is going to give rise to a specihc 
inflaton potential which we will derive below. 


Assuming flat FRW universe, and taking SvrG 
equation. 



= 1 yields to the following Friedmann 


( 2 . 2 ) 


where H is the Hubble parameter and a is the scale factor. Also, one can obtain conserva¬ 
tion equation as 

p + 3H{p + p) = 0. (2.3) 


Using the relations (2.1) and (2.2) in the equation (2.3) gives the following scale factor- 
dependence for energy density [45], 


/ B C 

\ 1 -|- A o^(i+-4)(i+«) J ’ 


(2.4) 


where C is an integration constant. 

It is clear from the equation (2.1) that during the early universe p ~ Ap. And for the late 
time, one can see from the equation (2.4) that energy density approaches to a constant 
P (pq:^) since a S> 1. If we put this value for p in the equation (2.1) for pressure; 


p = A{ 


B 


1 + A 


) i+“ — B( 


B 


1 + A' 


l+a = — ( 


B 


l + A 


)i+“ = -P, 


(2.5) 


The above equation tells us that the modihed Chaplygin gas behaves like cosmological 
constant during the later stage of expansion independent of parameters A, B or a. 

Let us assume that the quintessence scalar held 4> is minimally coupled to gravity [46], 
and it is the scalar held which derives inhation. We would like this scalar held to obey the 
equation of state of the modihed Chaplygin gas. One can achieve this by assuming that 
energy density and pressure are given by the following relations respectively, 

p = + 

P=l^^- ( 2 - 6 ) 
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where V{4>) is the potential for the scalar field. Using equations (2.6) and (2.1) one can 
get (j) to have the following form 




(2.7) 


Now, using the equations (2.2), (2.4) and (2.7) we can obtain. 


dcj) 

da 


\/3(l + A) 


Ua /l + 


Cil+A) 


Q,3(l+A)(l+a) 


( 2 . 8 ) 


The equation (2.8) has the following solution, 


^ 2 ^ 3(1 + ^) 

4> = — -Trr- -ttt:: -r cosh 

3(1 + ^4) (1 + a) 


^Y^l + + (j)f, (2.9) 


where 0/ is final value (when a S> 1) of the scalar field. Conversely, we can express the 
scale factor in terms of the scalar field, 


a = 


B 


-(—1 + cosh^ 9) 


where we defined. 


{1 + A)C' 

Q ^ 3(1 + A){1 + a) 


3(l+A)(l + a) 


-</>/)■ 


( 2 . 10 ) 


( 2 . 11 ) 


2^3{1 + A) 

Having known the dependence of the the scale factor on the scalar field we can obtain the 
expressions for the scalar potential and the Hubble expansion parameter as a function of 
the scalar field. Using the relations (2.4), (2.6) and (2.10) we can find. 


U = Uo 


^4 - 1 , o ■ 

1---cosh^ 6 

A + 1 


2a 

(cosh0) i+“ 


( 2 . 12 ) 


where Vq = ^^^(j^) is a constant that depends on A, B and a. Then, using the 
relations (2.2), (2.4) and (2.10) one can obtain. 


H = Ho(cosh6')U4^ , 


(2.13) 


where the constant Hq has a particular dependence on A, B and a of the following form, 

, p 1 

Ho = ! which we define for convenience. In the next section we will use above 

result to obtain the tensor to scalar ratio and fix the free parameters using Planck 2015 
+ BICEP2 data. In section 4, having fixed the free parameters of the model, we will 
analyze the behavior of cosmological parameters such as Hubble expansion parameter and 
the deceleration parameters. 
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3 The tensor to scalar ratio 


The first constraint that we will make use of is the tensor to scalar ratio. Assuming slow- 
roll inflation r, the tensor to scalar ratio which measures the amount of the stochastic 
gravitational wave production during inflation, is given by, 


r ~ IQeh- 


(3.1) 


Here en is one of the slow-roll parameters which is defined as, 

en = 2 


dH\‘^ 

1 


Using the equations (2.13) and (3.2) we get the following expression, 


e/f = -(1 + A) tanh^ 6. 


(3.2) 


(3.3) 


Therefore the tensor to scalar ratio can be expressed in terms of the scalar field as 

r ~ 24(1 + A) tanh^ 0. (3.4) 


On the other hand, the current upper bound on r by BKP (BICEP Keck Planck [47]) is 
r < 0.09. The equation (3.4) tells us that, in the early universe 6 —>• —oo then tanh^(0) —)■ 1 
and we have r —)• 24(1 + A). It means that the asymptotic behavior of the tensor to scalar 
ratio is 24(1 +A), hence, the only free parameter that determines the tensor to scalar ratio 
is A; which will be fixed using observational data of BKP as. 


- 1 < A < -0.996. (3.5) 

In the case of upper limit of the BKP joint analysis one gets A = —0.996. This is an 
expected result as we discussed in the introduction section and is different from the pre¬ 
vious results. For example, authors of the Ref. [48] suggested 0 < A < 1.35, also there 
are several works which used A > 0. Although there are some works which found small 
negative values for the parameter A, at least with our knowledge, the value of A ~ —0.99 
is a new result. This is due to enforcing the scalar held to be the cause of the early time 
and late time acceleration. In Fig. 1 we can see the behavior of r by various values of a. 
Assuming (p an increasing function of time, the limit of r < 0.09 is achieved for A ~ —0.99. 
This proves our expectation of the asymptotic behavior of r, that it is independent of all 
parameters other than A. 


Let us look at another important observable quantity, the scalar spectral index Ug 
which is given by. 


ns ~ 1 - 4eH + 25 h , 


(3.6) 


where. 




2HeHj ’ 


(3.7) 
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Figure 1. The tensor to scalar ratio in terms of the scalar field for i? = 3, C = 1 and (j)f = 0. 
(a) a = 0.5, A = —0.996 (dotted blue), A = —0.997 (dashed green), A = —0.998 (solid red), 
A = —0.999 (dash dotted yellow), (b) A = —0.996, a = 0.1 (dotted blue), a = 0.5 (solid red), 
a = 1 (dash dotted yellow). 
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Figure 2. The scalar spectral index in terms of the scalar field for a = 0.5, ( 7=1 and cj)f = 0. 
A = —0.996 (dotted blue), A = —0.997 (solid red), A = —0.998 (dashed green), A = —0.999 (dash 
dotted yellow). 

is the other slow-roll parameter. The current bound on Ug obtained by Planck 2015 TT + 
low P analysis [47] is n<j = 0.968 ± 0.006. In Fig. 2 we draw Ug in terms of the scalar field 
and see that, at the early universe (p —oo, the observed upper bound is satisfied. 

Having a negative but bigger than -1 value for A would constrain H to be a positive 
parameter. This can be seen from the following; equation (2.13) shows us that since 
Hq = only a positive value of B gives a positive real Hubble parameter. 
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The usual choice for the value of parameter a is between zero and unity. As an example, 
which is one among the many, in Ref. [41] the case of positive a = 1/2 is considered and the 
model is called the reduced modified Chaplygin gas model. In late time limit one could get 
the pressure contribution of Chaplygin gas to dominate the total pressure only if pressure 
is inversely proportional to density, which is a decreasing function of time. Despite this 
fact, in Ref. [39] some cases of negative a were also discussed. The sign and the value of 
a will be constrained using the current value of the Hubble parameter. In the next section 
we study cosmological parameters using obtained parameter values. 

4 Cosmological parameters 

Using equation (2.9) we can obtain the behavior of the scalar field with the scale factor 
as illustrated in Fig. 3. We can see that 4> is an increasing function of the scale factor. 
With the choice of 4>f = 0 the value of the scalar filed goes from negative infinity to zero. 
It is clear from Fig. 3 that the evolution of cj) at the initial states are faster that the later 
states. This might be a troublesome situation, since the early universe inflation is achieved 
by having the scalar held slowly roll from the inhaton potential, in turn to be small 
compared to scalar potential. Here we can see from the Fig. 4 that the magnitude of V{(j)) 
is so big compared to that, the doesn’t help to catch the value of the inhaton 
potential during the early universe for a long period of time. 



a 


Figure 3. The scalar field in terms of the scale factor for A = —0.996, a = 0.5, (7 = 1 and <// = 0. 
B = 1 (dotted blue), B = 2 (solid cyan)i? = 5 (solid red), B = 8 (dashed green), B = 10 (dash 
dotted yellow). 

Then, using equation (2.13) we can investigate the evolution of the Hubble parameter. 
One can see that the Hubble parameter is a decreasing function of scalar field and yields 
to a constant at the late times. Using the current value of Hubble parameter {H ~ 
67kms~^Mpc~^ ~ 10“®^) one can also fix B and a. The asymptotic value of Hubble 
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Figure 4. cfp (solid red) and V{4)) (dotted blue) in terms of the scalar field for A = —0.996, 
a = 0.5, B = 3, C = 1 and = 0. 

parameter is found in section 2 to be • To make this number close to 10“®^ 

(Hubble parameter’s current value in natural units) one has two choices: 

• a > 0, H ~ for Chaplygin gas. 


• a < 0, B ^ 0{1) for barotropic fluid. 


However they don’t have any effect on r and behavior of other cosmological parameters. 
Having H in terms of 4> one can obtain the deceleration parameter q as, 

~b" 


H 

g- -1- ^ - -1- 


•^(1 + ^)p — 


dH 




(4.1) 


So, we represent behavior of q in terms of the scale factor by Fig. 5 in exact agreement 
with observations which tells us that the final value of the deceleration parameter is about 
-1. We can see that the deceleration parameter initially begins with a negative value and 
approaches to -1 towards the end (late time). Above results suggest that our model can 
explain both the early and late time acceleration with a single inflaton potential V{(p). The 
problem is ending it during the intermediate stages. One can imagine a scenario where (p 
decays into other scalar particles at low energies and give rise to reheating but eventually 
dominate the universe due to its form (constant energy density). 

In the next section we include bulk viscosity and discuss about consequences. 


5 Effect of varying bulk viscosity 

In this section we assume that the space-time is filled with modified Chaplygin gas having 
a bulk viscosity. In that case, the total pressure changes as the following: 

Peff=P + '^- (5.1) 
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Figure 5. The deceleration parameter in terms of the scale factor for A = —0.996, B = 3, C = 1 
and (/)/ = 0. a = 0.1 (dotted blue), a = 0.5 (solid red), a = 1 (dash dotted yellow). 

Here the bulk viscous pressure H is represented by the Eckarts expression which is propor¬ 
tional to the Hubble parameter with proportionality factor identified as the bulk viscosity 
coefficient, 

H = -3^i7. (5.2) 

The bulk viscous coefficient is assumed to have a power law dependence on the energy 
density, 

C = ^op", (5.3) 

where and n are positive constants. Possible solutions were discussed for this proposed 
model [49]. The simplest solution obtained for the special case of n = 1/2 is considered 
here to obtain the effect of bulk viscosity on observed quantities such as and r. In the 
case of re = 1/2 one can obtain [49] the following expression for the energy density; 

B C 

1 -{- A — v^^o)(i+“) 

Repeating the similar calculations of the previous section shows us that the presence of 
viscosity is not suitable for modified Chaplygin gas model and this is in accordance with 
the suggestion that the modified Chaplygin gas behaving as viscous generalized Chaplygin 
gas. We find that value of the re* increases due to bulk viscosity, while the tensor to scalar 
ratio decreases by viscous coefficient. And the problem of getting the initial and final 
inflation consistently by this modification still remains. 

6 Conclusion and discussion 

In this paper we presented a model in which a single scalar field inflaton is responsible 
for both the primordial and late time acceleration. We forced this scalar field to obey the 
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modified Chaplygin gas EoS. By doing so, we constructed the inflaton potential bottom-up. 
In a recent work [50] a similar approach was used to explain both periods using a k-essence 
non-minimally coupled scalar field. Unlike giving the form of the potential initially; we 
wanted the Hubble parameter dictate the form of the inflaton potential by assuming Chap¬ 
lygin gas EoS via Friedmann equations. Recently [39], using Planck-I-BICEP2 data, it is 
found that generalized Chaplygin gas is not a suitable model of inflation in Einstein gravity, 
so one may need generalized Chaplygin gas inflationary model in RS type five dimensional 
brane world scenario. In this paper, we investigated the same problem for the modified 
Chaplygin gas model with and without viscosity in the light of more recent data [47]. In 
the modified Chaplygin gas model there are three fundamental parameters: A, B and a. 
We obtained that the most relevant free parameter to primordial inflation, A should lie 
between —1 < A < —0.996. The values for B and a are determined by the current value of 
Hubble parameter. The conclusion is that a is either a positive 0(1) number or a negative 
number which is very close to -1. The former case gives an unnaturally small value for 
B ~ or so. The latter case works with B ~ 0(1) which would certainly be pre¬ 

ferred. Therefore, one can say that late time acceleration caused by the same inflation field, 
prefers barotropic fluid EoS over Chaplygin gas EoS. We conclude that the presence of bulk 
viscosity is not necessary in modified Chaplygin gas and may yield to inconsistencies with 
observational data. 

The only shortcoming of our model is a mechanism to end inflation between two accelerated 
stages. This model which basically behaves quite similar to cosmological constant will ob¬ 
viously not allow production of large scale structure. This feature can be analyzed formally 
by solving perturbed Einstein’s equations but Fig. 5 leaves no room for that effort. Ending 
inflation can be achieved by some coupling to a spectator scalar field and therefore a decay 
to that particle would soon give rise to reheating. But this ad-hoc procedure should be 
developed more formally. 

For future work, it can be considered to use extended models including other free param¬ 
eters [22]. One can actually imagine a density term, with a postive constant in front, in 
the equation of state. If the power of that density term is negative than it would dominate 
the pressure between the two accelerated stages. One could solve the Friedmann equations 
to construct Inflaton potential for this case as well. Cosmological perturbation analysis 
would help to constrain additional free parameters that will enter with the new term. We 
can also investigate k-essence generalization and even non-minimal coupling between scalar 
field and gravity. All of these attempts will complicate the problem but nevertheless will 
give us more space to explain observable phenomena. 
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